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LAMINAR  BOUNDARY  LAYER  ON  CONE  IN  SUPERSONIC  FLOW  AT  LARGE  ANGLE 

OF  ATTACK  1 

Bv  Franklin  K.  Moork 


SUMMARY 

The  laminar  boundary-layer  flair  about  a  circular  cone  at 
large  angles  of  attack  to  a  supersonic  stream  has  been  analyzed 
in  the  plane  of  symmetry  by  a  method  applicable  in  general  to 
the  flow  about  conical  bodies. 

.1/  the  bottom  of  the  cone,  reincity  profiles  were  obtained  show¬ 
ing  the  expected  tendency  of  the  boundary  layer  to  become  thinner 
on  the  under  side  aj  the  cone  as  the  angle  of  attack  is  increased. 

At  the  top  of  the  cone,  the  analysis  failed  to  yield  unique 
solutions,  except  for  small  angles  of  attack.  Beyond  a  certain 
critical  angle  of  attack,  boundary-layer  flow  does  not  exist  in 
the  plane  of  symmetry,  thus  indicating  separation.  This 
critical  angle  is  presented  as  a  function  of  Mach  number  ami 
cone  vertex  angle. 

INTRODUCTION 

The  supersonic  aerodynamics  of  pointed  bodies  has  con¬ 
siderable  current  interest  in  connection  with  the  design  of 
aircraft  and  missile  fuselages.  An  important  feature  of  the 
flow  about  such  bodies  is  the  behavior  of  the  boundary  layer 
and,  in  particular,  the  flow  separation  which  may  occur  along 
the  low-pressure  side  of  the  body  due  to  angle  of  attack. 
The  present  report  will  consider  the  development  of  the  lami¬ 
nar  boundary  layer  on  the  surface  of  a  right  circular  cone  at 
an  angle  of  attack  to  a  supersonic  stream  (see  fig.  1).  'Pile 
conical  configuration  may  be  considered  an  idealization  of 
the  nose  portion  of  a  supersonic  aircraft  fuselage. 

Outside  a  thin  boundary  layer  on  a  cone,  the  nonviscous 
supersonic  flow  (upon  which  the  boundary  layer  itself  de¬ 
pends)  is  ‘‘conical”  in  the  sense  that  physical  quantities 


(such  ns  velocity  and  pressure)  are  constant  along  any  ray 
proceeding  from  the  cone  apex.  The  description  of  this 
outer  flow,  contained  in  references  1  to  A  and  elucidated  in 
reference  4,  is  considered  adequate  for  the  purposes  of  this 
report,  but  subject  to  restrictions  which  will  be  discussed 
subsequently. 

In  figure  2  is  shown  qualitatively  the  circumferential 
pressure  distribution  on  the  cone  surface  predicted  for  var¬ 
ious  angles  of  attack  (see  ref.  4).  These  pressure  distribu¬ 
tions  depend  only  on  the  character  of  the  nonviscous  flow 
beyond  the  boundary  layer,  on  the  assumption  that  the 
boundary  layer  is  extremely  thin.  When  the  angle  of  attack 
is  very  small,  the  pressure  decreases  inonotonically  from 
the  bottom  of  the  cone  around  to  the  top.  For  larger 
angles  of  attack  there  appears  a  region  near  the  top  of  the 
cone  wherein  the  pressure  gradient  reverses  and  the  pressure 
increases  toward  the  top.  As  the  angle  of  attack  is  furf lie- 
increased,  this  region  becomes  greater  in  extent. 


F«;i  kk  2.  -Pressure  distriluitiua  around  cone  fur  various  range-.  of 
angle  of  attack. 


As  a  consequence  of  the  conical  nature  of  the  nonviscous 
flow,  it  is  shown  in  references  ">  and  <i  that  the  laminar 
boundary  layer  lias  parabolic  similarity  along  generators  of 
the  cone;  that  is,  velocity,  pressure,  and  density  inside  the 
boundary  layer  are  constant  along  any  parabola  (see  fig.  1) 
drawn  in  any  one  meridional  plane,  (plant1  passing  through 
the  body  axis).  Of  course,  circumferential  variation  of  these 
quantities  is  to  lie  expected  when  the  cone  is  at  angle  of 
attack. 

In  reference  7  the  effect  of  angle  of  attack  on  flic  laminar 
boundary  layer  is  analyzed,  in  the  limit  of  very  small  angle 
of  attack,  with  the  result  that  the  boundary  layer  tends  to  be 
thicker  on  the  lop  of  the  cmie  than  on  the  bottom  (fig.  if(a)). 


1  Supersedes  N  AC  A  TN  2844,  “laminar  Boundary  Layer  on  Cone  In  Supersonic  Flow  at  Large  Angle  of  Attack"  by  Franklin  K.  Moore,  1952. 
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Full  he  3.  -Cross-sectional  views  of  boundary  layer  on  com*  ai  various 
angles  of  atlaek. 

This  is  to  he  expected  since  the  Hu  id  near  t  lie  base  of  the 
boundary  layer  has  low  inertia,  is  therefore  inclined  to 
follow  the  direction  of  the  circumferential  pressure  gradient 
more  closely  than  is  the  outer  How,  and  thus  tends  to  drain 
away  from  beneath  the  cone  and  accumulate  near  the  top. 
Xo  separation  is  encountered  because,  for  small  angle  of 
attack,  the  pressure  gradient  is  always  favorable  (Hg.  2(a)). 

For  larger  angles  of  attack,  when  the  pressure  gradient 
reverses  direction  near  the  top  of  the  cone  (fig.  2(b)),  experi¬ 
ment  indicates  the  formation  of  boundary-layer  “lobes" 


/ 

/  Total -head  lube 


(fig.  3(b)).  When  the  angle  of  attack  is  further  increased, 
the  lobe  put  tern  finally  breaks  away  from  the  body  t<>  form 
a  vortex  street  dig.  3(c)). 

Recently  (ref.  X).  at  the  Lewis  laboratory  a  brief  experi¬ 
ment  was  carried  out  in  which  a  total-head  probe  was  placed 
near  the  surface  at  the  top  of  a  cone  and  pointed  toward 
the  cone  apex  dig.  4(a)).  The  cone  was  mounted  in  a 
supersonic  wind  tunnel,  and  the  probe  was  used  to  measure 
the  total  head  ill  the  boundary  layer  at  a  fixed  height  above 
the  surface  as  the  angle  of  attack  was  varied  b\  rotating 
the  cone  in  the  meridional  plane  containing  the  probe. 
Figure  4(b)  shows  the  result  of  this  test.  The  decrease  in 
indicated  total  head  as  the  angle  of  attack  was  increased 
from  negative  to  positive  values  may  be  interpreted  to 
mean  an  increase  in  boundary-layer  thickness  at  the  top 
of  the  cone.  Beyond  a  certain  angle  of  attack,  this  tend¬ 
ency  reverses,  and  the  boundary  layer  apparently  becomes 
thinner  as  the  angle  of  attack  is  increased.  This  is  a  pos¬ 
sible  indication  of  the  tendency  to  form  lobes,  as  illustrated 
in  figure  3(b). 

In  the  present  report  of  research  conducted  during  the 
summer  of  H)”>2  at  the  \A('A  Lewis  laboratory,  the  laminar 
boundary  layer  in  the  meridional  plane  of  symmetry  of  the 
flow  is  analyzed  for  large  angles  of  attack  in  order  to  provide 
velocity  profiles  on  the  bottom  of  the  cone  and  to  provide 
a  certain  degree  of  insight  into  the  question  of  separation 
on  the  to]). 

SOLUTION  OF  BOUNDARY-LAYER  EQUATIONS  IN  PLANE  OF 
SYMMETRY 

BOUNDARY-LAYER  EQUATIONS  IN  FLANK  OF  SYMMETRY 

In  reference  7,  it  is  shown  that  the  dimensionless  laminar¬ 
houndary -layer  equations  for  supersonic  How  over  a  circular 
cone  are 


(t>)  Variation  of  total  head  with  angle  of  attack. 

Fie.cRK  4. — Variation  with  angle  of  attack  of  total  head  measured  in 
boundary  layer  of  c one.  M,  3.095;  0,  7.5°. 
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Equations  (la)  and  (lb)  are  momentum  equations,  equation 
(le)  is  an  energy  balance,  and  equation  (Id)  is  the  equation 
of  state.  A  complete  list  of  symbols  is  provided  in  appendix 
A.  The  functions  and  ,?(X,<£)  are  related  to  the  two- 

eomponent  vector  potential  discussed  in  reference  0  and  are 
defined  according  to  the  relations 


«  =  /x 

U'=.?X 


(2) 
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in  u  manner  such  as  to  satisfy  the  continuity  equation 
identically. 

The  coordinate  X  has  been  formed  as  follows: 


velocity  components  should  take  on  the  corresponding 
nonviseous  values 

./xt00 , vp)  ■  (tin) 


■ ’[(£,)  "1 /"'»} 


(3) 


Tlie  coordinate  p  is  the  angle  between  the  vertical  plane  of 
symmetry  of  the  flow  and  any  meridional  plane  of  the  body 
(fig.  1).  Equations  (1)  and  (If)  imply  that  parabolic  simi¬ 
larity  of  the  Blasius  type  exists  in  meridional  planes.  As 
pointed  out  in  references  5  and  fi.  this  conclusion  of  parabolic 
similarity  applies  for  the  boundary  layer  on  any  smooth 
conically  symmetric  body  in  supersonic  flow  (for  example,  a 
cone  of  elliptic  cross  section). 

In  reference  7,  all  quantities  are  made  dimensionless  by 
referring  them  to  the  properties  of  the  nonviseous  flow  at 
the  outer  edge  of  the  boundary  layer  when  the  cone  is  at 
zero  angle  of  attack.  In  the  present  report  it  will  be  con¬ 
venient  to  use  a  different  reference  condition  (subscript  r) 
which  will  be  defined  subsequently.  The  following  quan¬ 
tities  on  the  left  are  to  be  identified  with  the  dimensionless 
groups  on  the  right: 


>/,»'  -  >/•>/,,  »'  iir 

P  "  p‘p, 

T  -  2c ► 

P^P  PrU  r4 

x.y  -  ppi,i  ppipjjCn, " 


(4) 


where  the  constant  ('  arises  from  the  assumption  of  the 
temperature-viscosity  relation  of  Chapman  and  Rubesin 
(ref.  9): 


T 

Tr 


(5a ) 


with  C  being  defined  as  follows,  in  order  to  match  equation 
( 5a )  to  the  Sutherland  formula  at  the  cone  surface  (denoted 
by  subscript  w): 


fi—(  Tw\l  1  4-S/T, 

“V  Tj  T.lTr+S/Tr 


(5b) 


The  quantity  »S  may  be  taken  as  equal  to  (2 lti°  H)2rp/u/. 

The  following  additional  physical  assumptions  are  em¬ 
bodied  in  equations  (1): 

(a)  A  thin  boundary  layer  across  which  the  static  pressure 
is  constant 

(h)  Frandtl  number  of  1  and  constant  ratio  of  specific 
heats  7 

(c)  No  heat  transfer  through  the  surface 

From  equation  flc),  since  the  case  of  Fraud ti  number  of  1 
and  no  heat  transfer  is  considered.  T„  in  equation  (5b)  may 
he  taken  equal  to  the  dimensionless  stream  stagnation 
temperature. 

The  boundary  conditions  on  the  functions /(X,^)  and  ;i(\.p) 
are:  At  the  outer  edge  of  the  boundary  layer,  the  u  and  w 


.gx(°° ,  <^)  u\(p)  (Obi 

At  the  cone  surface,  the  u  and  >r  velocities  should  vanish 


/x(0,iP)^£x(<W>^-d  «>c) 


and  the  normal  velocity  r  should  vanish.  It  is  shown  in 
reference  fi  that  this  last  requirement  is  met  if 

/( ;g(0,v5l  —0  (Oil) 

Equations  (1 ).  involving  two  independent  variables,  would 
be  quite  difficult  to  solve  in  general.  However,  a  certain 
amount  of  information  can  be  obtained  by  restricting  con¬ 
sideration  to  the  plane  of  symmetry,  thus  yielding  a  tractable 
set  of  ordinary  equations  involving  X  as  the  only  independent 
variable. 

In  the  plane  of  symmetry  (^=0.ir),  n’=gx~  (I.  Because 
n  —  fk  is  even  about  the  plane  of  symmetry  and  may  he 
expected  to  be  regular  there.  /x#,  =  0.  The  pressure  and 
the  density  arc  also  even,  and  therefore  p'(p\  vanishes  at 
p— 0,7t.  Thus,  in  the  plane  of  symmetry,  equation  flat 
reduces  to  the  following  equation: 

^  /T  ~q  j  /,u  *  2/xxx  0  (7a) 


Every  term  in  equation  (lb)  vanishes  at  the  plane  of  sym¬ 
metry;  and,  therefore,  in  order  to  obtain  a  meaningful 
equation,  it  is  necessary  first  to  differentiate  equation  (lb) 
with  respect  to  p  and  then  drop  terms  which  vanish  at 
ip—  0.7T.  This  procedure  yields  the  following  result: 


(2  \  2  .,2 

.(  “b  g*  J  .g xx, fi  hq  0»x<p)  ^  5xip./x 


2  //'(*>) 
30  p 


-F  — gxxxv-  ~  1 1  (7b) 


Equation  (le)  becomes 


(7c) 


Equations  (7)  may  be  considered  a  set  of  ordinary  differ¬ 
ential  equations  for  the  functions  /(X.0)  and  g*(X.()),  or 
AX.ir)  and  depending  on  whether  the  solution  is 

required  at  the  bottom  or  the  top  of  the  cone.  According 
to  equations  (2),  the  result /tX.O  or  ir)  may  be  differentiated 
with  respect  to  X  to  give  the  profile  of  meridional  velocity  u 
in  the  plane  of  symmetry.  The  form  assumed  by  the  cir¬ 
cumferential  velocity  profile  «•  as  <p ■->()  or  ir  is  given  by 
g«,ji(X,0  or  ir)  in  the  sense  that,  at  a  small  angular  distance  dp 
away  from  the  plane  of  symmetry.  w^gf\  dp. 

The  boundary  conditions  (eqs.  (6))  become,  in  the  plane 
of  symmetry, 


/*(  00  ,0  or  ir)  —  u, (0  or  ir) 
gr\(  ®  ,0  or  ir)  =«',„(<)  or  ir) 

/„((),()  or  ir)  =  g^(0,0  or  ir)  =/((),()  or  ir)  =,§,(().()  or  ir)  =  0. 


(8) 


In  view  of  the  first  of  equations  (8)  and  of  equation  (7c).  it  is 
convenient  to  specify  the  reference  condition  (subscript  /•) 


4 


REPORT  1132  NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


lo  lx*  that  existin';  at  the  outer  edge  of  the  boundary  layer, 
for  the  particular  untile  of  attack  under  consideration, 
evaluated  at  either  (£•=()  or  <p=v,  dependin';  on  whether  the 
analysis  pertains  to  the  bottom  or  top  of  the  cone. 

Because  the  pressure  is  assumed  constant  across  the 
boundary  layer,  equations  (Id)  and  (7c)  and  the  assumption 
of  constant,  pressure  across  the  boundary  layer  (p=p,) 
provide  that 

X-  L  =  \  +  }r  [1-f^l  (9) 

From  conventions  (4).  f\  =  ( ~  ^  l/,'.  For  convenience, 

the  followin';  definitions  are  made: 

?,(X.O  or  =  k+{\)  (10a) 

£  =  30  w\(0  or  tt)  (10b) 


Equation  (9)  and  definitions  (10)  are  introduced  into  equa¬ 
tions  (7a)  and  (7b),  and  a  value  of  p"( 0  or  r)  is  assigned  con¬ 
sistent  with  the  nonviscous  equations  at  the  outer  edge  of 
the  boundary  layer  (/,"(()  or  ir)  may  conveniently  he  ob¬ 
tained  from  equations  (lb)  and  (ti)  by  setting  gvs\—gr\w^:0 
when  X=  =° ). 

The  following  pair  of  simultaneous  ordinary  differential 
equations  then  results: 

(/  +  £*)/"+2/"'=0  (11a) 

( £ + !)  j 1  +  yl-t  ( i — (/')21 1 = o  (in,) 

and  boundary  conditions  (8)  become 

/'( 00  )  =  t£'(  00  )  =  1  (12a) 

r(l))=t£'(0)  =  ()  (12b) 

,/(9) (C)=0  (12c) 

Two  parameters  appear:  k,  which  depends  essentially  on 
angle  of  attack,  and  1\,  which  is  essentially  dependent  on 
Mach  number.  If  the  angle  of  attack  (and  hence  k)  is  zero 
or  nearly  zero,  equations  (11)  become  precisely  those  con¬ 
sidered  in  reference  7  and  may  be  solved  quite  readily,  since 
equation  (lib)  becomes  linear  and  the  solution  of  equation 
(11a)  is  well  known  as  the  Blasius  function.  When  k  differs 
substantially  from  zero  (moderate  or  large  angle  of  attack), 
equations  (11)  are  both  nonlinear  and  the  solutions  are 
interdependent.  For  any  particular  case,  when  onlv  the 
stream  Mach  number,  cone  vertex  angle,  and  angle  of  attack 
are  specified,  the  parameters  k  and  Tx  must  be  obtained  bv 
recourse  to  a  theory  of  the  outer  nonviscous  flow. 

OUTER  NONVISCOUS  FLOW 

In  references  1  to  3,  the  results  of  a  theory  of  nonviscous 
supersonic  flow  about  circular  cones  at  angle  of  attack  are 


tabulated.  The  case  of  zero  angle  of  attack  (ref.  1 )  is 
solved  exactly  in  the  sense  that  no  assumption  of  small 
vertex  angle  is  made.  The  equations  are  then  expanded  in 
powers  of  angle  of  attack  with  tile  use  of  the  zero-angle-of- 
attack  solution  as  the  first  approximation.  Terms  linear  in 
angle  of  attack  are  presented  in  reference  2.  and  terms  pro¬ 
portional  to  the  square  of  angle  of  attack  are  presented  in 
reference  3.  Reference  4  clarifies  the  application  of  the 
theory  to  the  computation  of  flow  conditions  at  the  cone 
surface. 

There  are  two  objections  to  the  use  of  this  theory  in  the 
present  application: 

(1)  Xeglect  of  terms  in  the  expansion  beyond  that  in¬ 
volving  the  square  of  angle  of  attack  may  lead  to  an  insuffi¬ 
ciently  accurate  representation  of  the  flow  at  the  large  angles 
of  attack  which  are  of  interest.  Unfortunately,  no  non¬ 
viscous  theory  is  available  that  treats  the  effect  of  angle  of 
attack  with  greater  precision.  In  reference  4,  a  comparison 
of  the  theoretical  and  experimental  pressure  distributions  is 
presented  for  a  cone  of  semivertex  angle  of  10°,  at  a  Mach 
number  of  2,  and  of  an  angle  of  attack  of  12.2°.  The  agree¬ 
ment  shown  is  very  good,  especially  since  the  angle  of  attack 
is  sufficiently  large  that  the  pressure  distribution  is  of  the 
type  shown  in  figure  2(h). 

(2)  In  reference  10,  Ferri  points  out  t lint  the  method  of 
expansion  used  in  references  2  and  3  is  improper  near  the 
cone  surface  and  leads  to  an  erroneous  form  of  the  entropx 
distribution  around  the  cone.  Therefore,  the  theory  cannot 
be  applied  if  the  vertex  angle,  the  tingle  of  attack,  or  tin* 
Mach  number  is  so  large  t hat  the  flow  may  not  la*  considered 
essentially  isentropie.  In  reference  7,  an  argument  is  pre¬ 
sented  to  the  effect  that  in  the  limit  of  infinitesimal  angle  of 
attack  the  presence  of  a  boundary  layer  ensures  that  the 
error  in  entropy  distribution  is  of  no  consequence  even  for 
large  cone  vertex  angles.  That  argument  in  no  way  applies 
to  the  present  analysis  because  the  angles  of  attack  con¬ 
sidered  are  not  infinitesimal.  For  the  purposes  of  this 
report,  the  use  of  references  2  and  3  in  their  present  form  is 
justified  only  in  cases  for  which  isentropie  flow  may  he 
assumed. 

According  to  reference  4,  the  velocity  components  at  the 
cone  surface  are,  using  the  notation  of  references  1  to  3  for 
quantities  tabulated  therein, 
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The  pressure  and  density  are 
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The  barred  quantities  are  those  pertaining  to  the  ease  of 
zero  angle  of  attack.  From  equations  (10b),  (14),  and  (14), 
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!  The  plus  or  minus  sign  refers  to  <p  0  or  tt,  respectively 
From  conventions  (4 1 .  equations  (14).  (lot,  and  (Mil.  an< 

the  result  of  reference  1  (hat  1  //’  Tr,(l  —  Ti:)  1  .!/*’, 
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In  figure  5  are  shown  4  and  I  7',  as  functions  of  a  for  a  cone  j  1  “  — - 

of  semivertex  angle  of  7.5°  and  a  stream  Mach  number  of  | 

4.1.  From  the  tabulations  of  reference  2.  it  may  be  inferred 
that  under  these  conditions  the  isentropie  assumption  leads  | 
to  errors  of  less  than  l  percent  in  quantities  proportional  C 

to  the  angle  of  attack. 

SOLUTION  OF  EOLATIONS  AT  v-=>0 

!  0 

Equations  (11)  have  been  solved,  subject  to  boundary 
conditions  (12).  at  (p  =  0,  and  various  angles  of  attack  for  a  1 
nine  of  semivertex  angle  of  7.5°  and  a  stream  Mach  number  ! 
of  4.1  for  which  the  values  of  k  and  1  T,  are  given  in  figure  5. 

The  computations  were  carried  out  by  Dr.  Lynn  Albers  of  the 
Lewis  laboratory  and  are  described  in  appendix  11.  The 
resulting  boundary -layer  profiles  of  meridional  velocity  u  and 
gradient  of  circumferential  velocity  dw  dip  are  shown  in  figure  *  ~u 
ti.  'Plie  curves  for  a  -  0  are  obtained  from  reference  7.  The 
profiles  show  clearly  that,  as  the  angle  of  attack  is  increased, 
the  boundary  layer  becomes  thinner  on  the  bottom  of  the  _0 
cone,  and  the  shear  stress  at  the  wall  increases. 

Skin  friction. — The  meridional  and  circumferential  com¬ 
ponents  of  the.  viscous  shear  stress  at  the  bottom  of  the  cone 
surface  mav  be  written  in  coefficient  form  as  follows:  -1 
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where  the  quantities  on  the  right  arc  in  dimensional  form. 
Application  of  equations  (Id),  (2),  (4),  (4),  (5),  and  (10a) 
yields 


Angle  of  ottack,  a,  deg 

Fiocre  5.— Variation  of  parameters  k  and  1  7)  with  angle  of  attack. 
.1/.  3.1;  O,  7.5°. 
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where 

-P'u'x 

‘  Ml 

Variation  of  these  skin-friction  coefficients  with  angle  of 
attack  is  shown  in  figure  7  for  a  particular  case. 

Displacement  thickness. — In  reference  11  it  is  shown  that 
the  displacement  thickness  A  for  a  cone  at  angle  of  attack  is 
the  solution  of  the  equation 


where 


(20) 


(21) 


At  #>  =  0,  Wi  =  0  and  equation  (20)  may  be  solved  directly 


S*-!-  Jc&p 

~T+T 


(22) 
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(’sing  equations  (10)  and  applying  ('Hospital's  rule  to 
evaluate  at  ^=0'  the  limit  of  the  ratio  ic/mt  appearing  in 
equation  (21 )  yield 

R-/?’X’[,-#'+r."-'r,’]-'x  (2Sb> 

Therefore, 

Vrft  *-/-{,;*  0'-T)+li  -/■)[■  +  (1  -I  /'>]}  A 

(23c) 

Figure  8  shows  the  variation  of  displacement  thickness  with 
angle  of  attack  for  a  particular  ease  and  again  illustrates  the 
expected  progressive  shift  of  the  boundary  layer  from  bottom 
to  top  as  the  angle  of  attack  is  decreased. 

LIMITATIONS  OF  METHOD  AT  »>  =  . 

Except  for  quite  small  angles  of  attack,  equations  (11) 
cannot  be  solved  at  the  top  of  the  cone  (<p=^ir).  Over  part 
of  the  range  of  angle  of  attack,  the  solutions  are  indetermi¬ 
nate;  and.  beyond  a  certain  angle  of  attack,  the  solutions 
do  not  exist  at  all.  These  properties  of  equations  (11)  will 
he  demonstrated  and  discussed  in  the  following  paragraphs. 

Asymptotic  forms  of  equations.  The  difficulties  just 
mentioned  may  host  be  inferred  from  the  asymptotic  forms 
of  equations  (11)  At  large  X.  From  equation  (12a)  it  is  clear 
that,  for  large  X,  /  and  ^  may  he  written  as  follows: 

/=X-M\X) 

V'  — X-f'F(X) 

where  F'  (  ®  )  =  ().  Substitution  into  equations  (11) 

yields  the  asymptotic  forms  for  large1  X 

(1  +  k)\F"  4-  2F’"=  0  (24a) 

( 1  +  fr)X*"  +  24''"—  2  (je  +  3)  *'  =  l  F'  [ 1  +  tx  (*  +  3  )] 

(24b) 


Indeterminate  solutions.- -Consideration  will  now  he  given 
to  the  problem  of  obtaining  the  complementary  solution  of 
equation  (24b).  Defining  a  new  dependent  variable 


L±* 

G  =  e  *  4'' 


yields  the  equation 

*+i§+( 1  ^7x2]r;==0  (25) 

This  is  essentially  Weber’s  equation  (ref.  12,  paragraph 
16.5),  and  the  asymptotic  solutions  are 

..  *  F  S/O 


k  +  1/3 

V+i 


t 

Thus,  the  asymptotic  complementary  solutions  of  equation 
(24b)  are 


Because  it  is  required  that  4',(  00  1  (|.  -olutiiui  1 2  *  >  I  *  ■  i- 
rejected  when  k —  1  .1.  If  -  1  5.  both  solutions  mav  he 

retained,  and  an  additional  undetermined  constant  appear- 
When  /•  >()  (at  the  bottom  of  the  conei,  the  complete  solu¬ 
tion  of  (‘((nations  (II)  i‘\ists  and  is  unique,  and  solution 
(26b)  is  to  be  rejected  in  forming  the  asymptotic  solution. 
Therefore,  if  solution  (26b)  and  the  associated  constant  must 
be  retained  when  /■<  I  5.  it  is  dear  that  the  complete 
solution  for  —  1  cannot  bo  unique.  This  iudctcrmimicx 
has  been  verified  by  numerical  integration  of  equations  i  1  I  1, 
as  described  in  appendix  B. 

This  indeterminacy  arises  because  essential  information 
has  been  lost  by  specializing  the  equations  to  apply  only  in 
the  plane  of  symmetry.  When  the  equations  are  so  written, 
it  is  implied  that  the  boundary-layer  development  in  the 
plane  of  symmetry  is  affected  only  by  conditions  in  that 
plane.  The  lateral  region  of  influence  of  points  on  the  bod\ 
in  the  plane  of  symmetry  grows  parabolicalh  (the  shaded 
regions  in  sketches  (a)  and  (bn,  according  to  the  law  of 
molecular  diffusion,  when,  as  in  the  present  instance,  there 
is  no  pressure  gradient  in  the  stream  direction.  Fluid  enters 
the  boundary  layer  from  the  outer  stream.  If  the  fluid  then 
moves  laterally  (because  of  the  angle  of  attack)  out  of  tin* 
region  of  influence  of  the  plane  of  symmetry,  sis  shown  in 
sketch  (a),  the  flow  is  uniquely  determined  by  outer  stream 
conditions  in  the  plane  of  symmetry.  Clearly,  this  is  the 
ease  when  <p  -  v  and  a  is  small,  and  when  yr  0  and  «  has  any 
positive  value.  When  except  for  -mall  angles  of 

attack,  the  lateral  motion  of  the  fluid  is  inward  relative  to 
the  region  of  influence,  as  shown  in  sketch  ( 1  > ) .  This  fluid 
(hen  brings  into  the  region  of  influence  of  the  plane  of  sym¬ 
metry  information  concerning  boundary-layer  development 
as  it  proceeds  around  the  cone  from  the  bottom.  Conse¬ 
quently,  outer  stream  conditions  at  <p  t  may  not  uniquely 
determine  the  boundary -layer  characteristics  at  #  -k.  and 

indeterminate  solutions  of  equal  ions  (11)  may  be  anticipated. 
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Tin1  condition  k~  —  specifies  an  angle  of  at  lack  such  (hat, 
at  <p~  w.  the  outer  flow  streamlines  are  just  tangent  to  the 
parabolic  region  of  influence.  This  statement  may  Im¬ 
proved  as  follows:  From  definition  (lOh).  if  k  —  b,  then 
■M*i  'Hi-  —  0  2.  .Vear  <p-  w  (see  sketch  (c)),  »'i  -  H’i/t 

W’i  -  •  Thus,  when  k—~!i.  i/v'iq-  r'2.  From  this  last 

(-(Illation  and  the  geometrical  properties  of  parabolas,  it  may 
be  inferred  that  the  streamlines  are  parabolas  with  focus  at 
the  cone  apex.  Therefore,  the  situation  shown  in  sketch 
(b)  applies  if  and  equations  (II)  have  no  unique 

solution. 

It  is  noteworthy  that  the  mathematical  nature  of  the 
indeterminacy  so  far  described  is  quite  similar  to  that  which 
characterizes  solutions  of  tin-  Falkner-Skan  equation  for  (lie 
plane  boundary  layer  when  the  outer  flow  velocity  is  pro¬ 
portional  to  x~n  (x  being  tin1  streamwise  coordinate  and  in  a 
positive  number).  Hartree  (ref.  Id)  has  treated  this  prob¬ 
lem  in  detail  and  finds  an  “extra”  asymptotic  solution  similar 
to  equation  (2(>b),  decreasing  with  distance  away  from  tin- 
surface  according  to  a  power  law.  This  extra  solution  is 
rejected  by  Hartree  on  the  ground  that  an  unrealistically 
thick  boundary  layer  would  be  predicted  if  the  solution 
were  retained. 

A  thorough  study  of  the  asymptotic  nature  of  equations 
0  )  might  lead  to  a  similar  conclusion  in  tlx-  present  problem. 
However,  in  the  present  problem,  the  mathematical  difficulty 
goes  deeper;  when  £>  — b,  numerical  integration  failed  to 
give  unique  results,  even  though  solution  (2tib)  is  certainly 
to  be  rejected.  Thus,  the  condition  k  >  — is  not  sufficient, 
but  rather  only  necessary  for  uniqueness. 

Inspection  of  the  profiles  shown  in  figure  fi  fora  =  k  —  0  indi¬ 
cates  that,  for  0>£>  — !il,  the  streamlines  w  ithin  the  bound¬ 
ary  layer  may  be  expected  to  incline  more  sharply  toward 
the  plane  tp=ir  than  do  the  streamlines  near  the  edge  of  the 
boundary  layer  and  thus  may  bring  information  from  beneath 
the  cone  even  though  the  outer  ones  do  not.  Therefore,  the 
necessary  condition  for  umqueu-ss  would  be  (see  eqs.  (2) 
and  (10a)): 


Figure  0  indicates  that  perhaps  the  maximum  value  of  ip'lf 
is  to  be  found  at  X--0,  in  which  case  criterion  (27a)  would 
become 


.  i/"(0) 

'  3*"(0) 


(27b) 


Nonexistence. —The  solution  of  equation  (24a)  is 


the  boundary  conditions,  and.  therefore,  the  Fremiti  boi  darv  - 
layer  equations  fail  to  describe  the  llow  This  was  first 
pointed  out  by  Hayes  in  reference  .V  The  Frandtl  equations 
differ  from  the  exact  equations  essentially  in  that  a  thin 
boundary  layer  is  assumed.  Thus,  if  k<  1.  the  boundarv 
layer  cannot  be  regarded  as  thin.  It  may  be  noted  that 
oqunt;>>n  (22)  implies  that,  as  k  ■  — !,  the  displacement 
thickness  approaches  infinity. 

Any  boundary  layer  grows  by  the  entrainment  of  fluid  at 
its  outer  edge.  That  is,  Iluid  particles  acquire  vortieitv  hv 
entry  into  the  boundary  layer.  Tin-  reverse  process  cannot 
occur  Iluid  particles  cannot  leave  the  boundar\  laver.  tlm- 
losing  their  vortieitv.  In  the  case  under  consideration,  it 
will  be  shown  that  when  k-  --  I.  the  streamlines  at  the  outer 
edge  of  the  boundary  layer  would  proceed  outward  relative 
to  the  boundary  layer,  if  the  boundary  laver  wore  to  retain 
parabolic  similarity.  Because  such  a  situation  is  phvsieallv 
impossible,  the  Frandtl  equations  fail  to  yield  a  solution. 

From  reference  II.  the  normal  velocity  at  the  outer  edge 
of  the  boundary  layer  at  p  it  is 
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d± 

Or 
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vvlicre  h  is  a  somewhat  arbitrary  definition  of  the  outer  edge 


of  the  boundary  layer  and  (  1  )  is  obtained  from  analysis 


of  the  outer  nonviscous  flow.  The  equation  of  enntimiitv 
for  the  outer  llow.  evaluated  at  the  surface  of  the  cone  in  the 
plane  of  symmetry,  may  be  w  ritten 
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I 

ex 
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With  equations  I.2N)  and  (21))  combined,  the  (low  inclination 
at  the  outer  edge  of  the  boundary  layer  is 


-J.  A  f  i  +  1  \ 

dx  x  \  '  Du,  dp  ' 


(30) 


The  Frandtl  (thin)  boundary  layer  may  be  expected  to  exist 
only  if 

dli  ^  r  i 
dx'  n  i 


or,  with  equation  (30)  introduced,  if 
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(31) 


If  parabolic  similarity  is  assumed  {It  and  A  each  proportional 
to  ^7),  inequality  (31)  becomes 


1  die,  _3 

diii  dip  ^  2 


The  requirement  that  F"{  «>)=0  is  met  only  if  £>—  1.  If 
k<  —  1,  no  solution  of  equation  (24a)  exists  which  satisfies 


or,  from  equation  (10b), 

£>-1 
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SEPARATION 

Till*  critical  condition  Jt‘=  —  1  might  he  expected  to  he  of 
physical  as  well  as  analytical  significance  because  some  sort 
of  catastrophic  thickening  of  the  houndary  layer  is  implied 
In  particular,  this  critical  condition  may  reasonably  tie 
supposed  to  he  connected  with  the  phenomenon  of  laminar 
separation.  In  order  to  explore  this  possibility,  it  is  first 
necessary  to  describe  qualitatively  what  is  meant  by  separa¬ 
tion  in  three-dimensional  boundary-layer  flow.  Difficulty 
has  been  encountered  in  establishing  a  satisfactory  qualitative 
criterion  for  three-dimensional  separation  (see,  for  example, 
ref.  a).  Therefore,  in  the  subsequent  paragraphs,  the 
general  problem  of  three-dimensional  separation  will  be 
discussed,  and  then  the  particular  case  of  the  cone  at  angle 
of  attack  will  lie  considered. 

CKNBRAL  CONSIOKR ATIONS 

Plane  flow.  In  plane  flow,  separation  is  customarily 
identified  by  t lu*  appearance  of  reverse  flow  (sketch  («!)). 
In  order  to  generalize  this  concept  to  three-dimensional 
How.  it  is  necessary  to  consider  the  separated  region  as  a 
whole.  In  plane  flow,  the  separation  point  of  sketch  (d.) 
might  be  regarded  as  the  forward  boundary  of  a  vortex 
sheet  embedded,  or  encapsulated,  within  a  region  bounded 
by  the  body  and  a  stream  surface  meeting  the  body  (sketch 
(eii.  Sketch  (c;  show-  separation  followed  by  reattach¬ 
ment.  Of  course,  the  sort  of  separation  of  greatest  engineer¬ 
ing  importance  occurs  when  such  an  embedded  vortex  sheet- 
rolls  tip  to  form  a  large  concentrated  vortex,  or  is  shed  as  a 
vortex  street,  with  the  consequence  that  the  outer  flow  is 
greatly  disturbed  and  a  large  pressure  effect  (form  drag) 
occurs. 


If  a  boundary-layer  solution  of  the  type  shown  in  sketch 
(c'  were  obtained,  the  vortex  sheet  would  be  completely 
embedded  in  a  thin  boundary  layer,  and  would  presumably 
tend  to  remain  flat  against  the  body  (leaving  out  of  account 
possible  effects  of  laminar  instability).  However,  if  the 
theoretically  predicted  vortex  sheet  extends  downstream 
into  a  region  where  the  thin  boundary-layer  equations  do 
not  apply  (that  is,  where  the  solution  “blows  up"  predicting 
an  infinitely  thick  boundary  layer),  then  in  that  region  the 
vortex  sheet  would  not  be  constrained  to  lie  flat  and  the 
rolling-up  process  would  occur. 

The  foregoing  discussion  seems  to  provide  the  proper 
description  of  what  occurs  behind  a  bluff  body:  A  complete 
solution  of  the  thin  laminar  boundary -layer  equations  for 
the  flow  over  a  cylinder  would  probably  yield  a  streamline 


pattern  of  the  form  shown  in  sketch  (ft.  The  botmdari 
layer  would  lie  predicted  to  gain  mass  flow  by  entrainment 
as  it  proceeds  around  the  body,  until  it  reaches  the  vicinity 
of  the  rear  stagnat ion  point .  There,  the  mass  flow  contained 
in  the  boundary  layer  must  finally  leave  the  body  and  proceed 
downstream.  The  boundary  layer  therefore  cannot  remain 
theoretically  thin,  but  rather  must  approach  infinite  thick¬ 
ness  in  violation  of  the  I’randtl  assumptions.  In  this 
region,  then,  the.  aft  boundaries  of  the  pair  of  vortex  sheets 
are  free  to  roll  up  into  concentrated  vortices,  thus  distorting 
the  outer  (low  in  such  a  way  that  the  rolling-up  process 
engulfs  most  of  the  region  which  would  otherwise  lie  occupied 
by  a  flat  vortex  sheet.  The  leading  edge  of  the  sheet, 
however,  is  still  constrained  to  lie  flat  against  the  body. 

Three-dimensional  flow.  The  foregoing  description  of 
plane  separation  may  be  generalized  to  three-dimensional 
flow  ns  follows:  A  separated  region  on  a  three-dimensional 
body  consists  of  a  vortex  sheet  embedded  between  the  body 
surface  and  a  stream  surface  attached  to  the  body  in  a  closed 
curve,  as  shown  in  sketch  (g),  which  is  a  view  of  the  body 
from  above.  'Pile  arrows  indicate  possible  directions  of 
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resultant  shear  stress  at  the  surface  aiul  outside  the  separated 
region.  The  situation  shown  in  the  sketch  would  correspond 
to  separation  and  reattachment  in  plane  flow.  If  somewhere 
within  or  at  the  boundary  of  the  separated  region  the  theo¬ 
retical  boundary  -layer  solution  would  a  priori  be  expected 
to  blow  up.  then  the  vortex  system  within  the  separated 
region  is  free  to  roll  up  into  a  more  or  less  vigorous  system  of 
vortices. 

Thus,  thin-boundarv-laver  theory  can  be  used  to  obtain 
the  following  information  concerning  laminar  separation: 

(1)  The  solution  may  establish  the  existence  of  a  vortex 
sheet  which  is  embedded  in  a  flat  bubble  on  the  surface  and 
which  could  adhere  to  the  surface  ami  remain  part  of  a  thin 
boundary  layer,  provided  that  the  Prandtl  equations  are 
valid  everywhere  in  the  separated  region. 

(2)  The  solution  might  predict  the  boundary  layer  to  go 
to  infinite  thickness  somewhere  in  the  separated  region.  If 
this  happens,  then  the  separated  region  is  free  to  roll  up,  thus 
providing  a  vigorous  wake  (which,  of  course,  is  not  amenable 
to  boundary-layer  theory). 

For  flow  about  a  plane  body,  the  boundary -layer  solution 
is  not  needed  for  predicting  the  breakdown  of  tbe  I’rnndtl 
assumptions.  Physical  considerations  suffice  to  establish 
when*  (at  stagnation  point  of  outer  flow)  and  when  (always) 
the  breakdown  occurs.  In  three-dimensional  flow  this  is  not 
always  so  clear,  and  at  least  certain  features  of  the  solution 
are  required  to  be  known.  In  order  that  the  boundary-layer 
equations  be  applicable,  the  solution  must  be  such  that  the 
boundary  layer  entrains  fluid  (that  is,  flow  streamlines  enter, 
but  do  not  leave,  the  boundary  layer  at  its  outer  edge).  In 
a  ('artesian  system,  where  A(j\r)  is  the  outer  edge  of  the 
boundary  layer,  this  requirement  may  be  written  as 

(/,/(,+ J/'i/^XriV*  (32) 

Equations  (32)  and  ( 2S )  may  be  combined  with  the  equation 
of  continuity  in  the  form 


to  yield 

|(A-A)piU,]-i-dv  |(A— Alp,ir,|>0 

or.  in  vector  notation, 

div  [(A  —  A) p i </ , | > 0  (33) 

where  A  is  the  displacement  thickness,  and  <j\  is  the  velocity 
vector  in  the  outer  flow  evaluated  at  the  body  surface. 


In  many  cases,  circumstances  may  be  found  for  which 
inequality  (33)  cannot  be  satisfied.  For  example,  for  plane 
incompressible  flow  about  a  cylinder  and,  as  is  customary, 
with  i/  and  x  defined  parallel  to  the  surface,  inequality  (33) 
may  be  written 


tlx 


(h -A) 


A-A 


du, 

tlx 

Ui 


As  the  rear  stagnation  point  of  the  outer  flow  is  approached. 
it i  tends  to  zero,  while  —<///,  tlx  remains  finite  and  positive. 
Therefore,  since  A  — A  by  definition  must  be  greater  than 
zero,  //(A— A)  tlx  must  approach  infinity  in  clear  violation 
of  the  boundary-layer  assumptions. 

SEPARATION  ON  CONE  AT  ANUI.E  OF  ATTACK 

In  the  previous  discussion,  it  was  concluded  that  separation 
involving  a  strong  vortex  pattern  occurs  if  a  tentative 
boundary-’nver  solution  predicts  an  embedded  vortex  sheet 
coupled  with  a  local  breakdown  of  the  assumption  of  a  thin 
boundary  layer.  In  the  case  of  the  cone,  incquulit  v  (33  < 
may  be  used  to  predict  the  circumstances  under  which  the 
boundary  layer  may  not  be  regarded  as  thin:  When  tbe  fact 
of  parabolic  similarity  is  introduced  (A  —  A  proportional 
to  —\x),  inequality  (33)  becomes 


d 


(A  -A) 


-  »'  i 


* 


3 

■> 


6h. 


A— -A  if  i 


(34) 


This  inequality  indicates  infinite  boundary -layer  thick¬ 
ness  only  when  if,  -  0  and,  even  then,  only  if  iru.  is  negative 
(which  is  true  at  the  top  of  the  cone,  p  it)  and  larger  in 
magnitude  than  ?,9it This  is  true  only  for  angles  of  attack 
larger  than  that  for  which  A- —  1  (by  eq.  (lOb)t.  When 
the  angle  of  attack  is  smaller  than  this  critical  value,  the 
right  member  of  inequality  (34)  is  always  negative,  and 
d(A—  A)'cV  may  be  considered  to  vanish  by  symmetry 
at  tp ir  without  violating  inequality  (34). 

The  foregoing  result  may  be  explained  on  physical  grounds 
as  follows:  As  the  boundary  layer  proceeds  around  the  cone, 
it  entrains  fluid  which  it  then  conveys  toward  the  top 
(symmetrically,  from  both  sides  of  the  cone).  In  the  plane 
cylinder  case,  the  fluid  similarly  conveyed  must  finally 
erupt  from  the  boundary  layer  when  the  stagnation  point 
is  reached.  However,  on  the  cone,  the  boundary  layer  grows 
parabolically  along  generators;  and.  hence,  if  the  crossflow 
is  not  too  strong  (small  angle  of  attack),  the  fluid  brought 
to  the  top  may  simply  become  part  of  the  growing  boundary 
layer.  For  larger  angles  of  attack,  the  boundary  layer 
cannot  grow  at  a  rate  sufficient  to  absorb  the  additional 
fluid,  and  eruption  occurs  with  tin'  consequent  breakdown 
of  the  thin-boundary-layer  assumptions. 

Accordingly,  it  is  proposed  that  when  the  angle  of  attack 
is  less  than  that  for  which  A'  —I,  a  thin  boundary  layer 
may  cover  the  cone  (fig.  3(a)).  For  larger  angles  of  attack, 
any  vortex  sheet  present  will  roll  up  to  form  attached  lobes 
(fig.  3(b));  for  still  larger  angles  of  attack,  a  vortex  street 
is  produced. 

Thus,  when  A-<— 1  (angle  of  attack  greater  than  that  for 
which  A-— 1),  strong  viscous  cross  forces  (viscous  lift)  on  the 
cone  may  be  expected.  These  forces  are  discussed  by  Allen 
and  Perkins  in  reference  14.  Of  course,  a  weakness  of  the 
present  analysis  is  that  no  indication  is  given  as  to  the 
strength  of  the  rolled-up  vortex  system  because,  when  A-  —  —  1 , 
the  presence  of  an  embedded  vortex  sheet  over  the  top  part 
of  the  cone  has  not  been  established.  It  seems  likely  that 
such  a  vortex  sheet  does  exist  because  a  rather  strong  adverse 
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f—r. 

pressure  gradient  (fig.  2(h))  always  exists  when  k—  —  1.  In 
fact,  it  may  he  shown  (most  conveniently  by  evaluating  eq. 
(7a)  at  as)  that  an  adverse  pressure  gradient  exists 
when  /'<  —2/3. 

Equations  (13),  (14),  and  (1  Oh)  give  the  critical  angle  of 
attack  corresponding  to  k~  —  1.  This  angle  is  presented  as 
a  function  of  stream  Mach  number  and  vertex  angle.  Figure 
9  shows  the  results  of  such  a  calculation.  The  critical  angle 
of  attack  is  given  as  a  ratio  of  angle  of  attack  to  semivertex 
angle  for  convenience.  The  results  suggest  that,  in  general, 
separation  involving  lobes  occurs  later  (in  terms  of  relative 
angle  of  attack  a/B)  for  the  smaller  vertex  angles,  particu¬ 
larly  at  higher  Mach  numbers.  Figure  9  indicates  the  pos¬ 
sibility  of  rather  profound  qualitative  differences  in  the 
flows  at  high  Mach  number  about  cones  of  different  vertex 
angles. 

Tht>  foregoing  interpretation  of  the  critical  condition 
k~  —  1  is  supported  by  the  experimental  result  shown  in 
figure  4.  From  figure  5.  k—  —  1  when  a=t>.2°,  under  the 
conditions  of  the  test.  Figure  4(h)  shows  the  measured 
total  head  rising  as  the  angle  of  attack  is  increased  beyond 
(i.2°.  Possibly  this  effect  is  caused  by  the  induced  field  of 
the  symmetrical  pair  of  vortex  lobes  sweeping  away  the 
thick  boundary  layer  between,  thus  reestablishing  a  thin 
boundary  layer  at  the  top  of  the  cone. 

It  may  he  of  interest  to  note  that  if  similarity  also  holds 
for  the  turbulent  boundary  layer  on  a  cone,  and  the  similarity 
law  is  nearly  linear  (rather  than  parabolic  as  in  the  laminar 
case),  separation  would  first  appear  at  a  higher  angle  of 
attack  than  in  laminar  flow.  In  fact,  equations  (!0b)  and 
(31)  or  (33)  would  yield  the  criterion  k=  —  4/3. 

For  the  boundary  layer  produced  by  supersonic  flow  over 
any  smooth  conically  symmetric  body  in  supersonic  flow 
(such  as  a  cone  of  elliptic  cross  section),  inequality  (33)  and 
the  condition  of  parabolic  boundary -layer  similarity  may  be 
used  to  find  a  criterion  equivalent  to  £=  —  1  for  the  maximum 
angle  of  attack  consistent  with  a  thin  boundary  layer. 


CONCLUSIONS 

The  laminar  boundary-layer  flow  about  a  circular  cone  at 
large  angles  of  attack  to  a  supersonic  stream  has  been  ana¬ 
lyzed  in  the  plane  of  symmetry  with  the  following  results: 

1.  At  the  bottom  of  the  cone,  profiles  of  meridional  ve¬ 
locity  and  of  the  gradient  of  circumferential  velocity  were 
determined  and  showed  the  expected  tendency  of  the  bound¬ 
ary  layer  to  become  thinner  on  the  underside  of  the  cone  us 
the  angle  of  attack  is  increased. 

2.  At  the  top  of  the  cone,  except  for  very  small  angles  of 
attack,  the  analysis  (which  is  restricted  to  the  plane  of  sym¬ 
metry)  failed  for  the  following  reasons: 

(a)  For  angles  of  attack  greater  than  some  rather  small 
value,  the  boundary  layer  brings  information  from  beneath 
the  cone  into  the  vicinity  of  the  plane  of  symmetry  at  the  top. 
Therefore,  the  analysis,  which  deals  only  with  the  plane  of 
symmetry,  yielded  indeterminate  solutions. 

(b)  For  angles  of  attack  greater  than  some  angle  (roughly 
of  the  order  of  the  cone  semivertex  angle),  no  boundary- 
layer  solution  is  possible.  The  characteristics  of  the  outer 
flow  and  the  known  parabolic  similarity  of  the  boundary 
layer  would  together  imply  that,  beyond  this  critical  angle, 
there  would  be  a  component  of  flow  leaving  the  boundary 
layer.  This  is  physically  impossible,  since  a  boundary  layer 
always  entrains  fluid.  Thus,  beyond  the  critical  angle  of 
attack,  no  solution  can  exist  for  equations  which  presume  a 
thin  boundary  layer. 

For  three-dimensional  flow  it  is  proposed  that  a  separated 
region  be  regarded  as  a  vortex  sheet  embedded  in  the  bound¬ 
ary  layer,  remaining  flat  against  the  body  if  the  assumption 
of  a  thin  boundary  layer  is  valid  throughout  the  region.  If, 
however,  the  boundary-layer  assumptions  break  down  any¬ 
where  in  the  separated  region,  it  is  inferred  that  the  vortex 
sheet  may  roll  up  to  form  strong  vortices  which  may  either 
remain  attached  or  be  shed  as  a  vortex  street. 

On  the  cone,  therefore,  the  critical  angle  of  attack  beyond 
which  no  boundary -layer  solution  is  possible  at  the  top  of 
the  cone  represents  the  maximum  angle  of  attack  for  which 
the  boundary  layer  is  everywhere  thin  or,  alternatively,  the 
minimum  angle  of  attack  for  which  major  disruption  of  the 
flow  may  be  expected  because  of  the  formation  of  strong 
vortex  lobes.  Beyond  this  angle  of  attack,  strong  viscous 
cross  forces  may  be  anticipated. 

A  similar  criterion  could  easily  be  obtained  for  the  bound¬ 
ary  layer  on  any  smooth  conically  symmetric  body  in  super¬ 
sonic  flow. 

The  assumption  of  a  suitable  similarity  law  suffices  to  es¬ 
tablish  a  similar  criterion  if  the  boundary  layer  on  a  conical 
body  is  turbulent. 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  September  15,  W52 


12 


REPORT  1  132 —  NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


APPENDIX  A 


SYMBOLS 


The  following  symbols  are  used  in  this  report: 

C  constant  appearing  in  temperature-viscosity 

relation  (eq.  (oa)) 

C,  component  of  skin-friction  coefficient  in 

/•direction 


o. 


F(\) 

KKv) 

h 

k 

M 

P 

<h 

nz 

Ri.  i\ 

T 

u 

V 

w 

x 

y 


component  of  skin-friction  coefficient  in 
^-direction 

specific  heat  at  constant  pressure 
function  appearing  in  asymptotic  representa¬ 
tion  of  /  (eqs.  (24)) 

function  related  to  meridional  velocity  u  by 
eq.  (2) 

function  related  to  circumferential  velocity 
w  by  eq.  (2) 

height  of  outer  edge  of  boundary  layer 
related  to  circumferential  gradient  of  cir¬ 
cumferential  velocity  in  plane  of  symmetry 
(eq.  (10b)) 

Mach  number 
static  pressure 

velocity  vector  at  outer  edge  of  boundary 
layer 

Reynolds  number,  pi«tx/gi 
Reynolds  numbers,  p,u iA/p,,  poh&i/iM, 
PiWi5»/mi(  respectively 
absolute  static  temperature 
meridional  component  of  velocity 
component  of  velocity  normal  to  surface 
circumferential  velocity  component 
coordinate  along  generators  of  cone 
coordinate  normal  to  surface 
angle  of  attack  (positive  as  shown  in  fig.  1) 


y 

A 

St 

0 

e 

x 

M 

P 

¥> 

*(x) 

*(X) 

Subscripts: 

max 

r 


1 


Superscripts: 


ratio  of  specific  heats 
displacement  surface  height 
mass-flow  defect  associated  with  meridional 
velocity  profile  (eq.  (21)) 
mass-flow  defect  associated  with  circum¬ 
ferential  profile  (eq.  (21)) 
semivertex  angle  of  cone 
sine  of  semi  vert  ex  angle'  of  coni' 
dimensionless  variable  (eq.  (2)) 
coefficient  of  viscosity 
density 

angular  coordinate  around  cone 
function  appearing  in  asymptotic  representa¬ 
tion  of  (eqs.  (24)) 

function  related  to  circumferential  velocity 
w  in  plane  of  symmetry  by  equation  (10a) 

maximum 

reference  condition,  non  viscous  flow  at 
surface,  at  <p  —  l)  or  <p  =  ir,  whichever  is 
appropriate 

evaluation  at  outer  edge  of  boundary  layer 
(alternatively,  non  viscous  flow  at  surface) 
Subscript  notation  for  partial  differentiation 
has  been  used 

Primes  denote  ordinary  differentiation  with 
respect  to  X  or  <p 

liar  over  quantity  indicates  evaluation  of 
nonviscous  flow  at  surface  when  cone  is 
at  zero  angle  of  attack 


APPENDIX  B 

NUMERICAL  SOLUTION  OF  DIFFERENTIAL  EQUATIONS 

By  Lynx  Albeks 


The  two  simultaneous  nonlinear  ordinary  differential 
equations  (Ua)  and  (lib)  together  with  boundary  conditions 
(12)  constitute  a  two-point  boundary-value  problem.  The 
method  of  numerical  solution  used  applies  directly  only  to 
problems  for  which  all  boundary  conditions  are  specified  at 
a  single  initial  point  (the  origin,  in  the  present  case).  Each 
numerical  integration  was  therefore  performed  starting  with 
boundary  conditions  (12a)  and  (12b)  and  a  tentative  speci¬ 
fication  of  /"(0)  and  ^"(0).  In  each  case,  such  integration 
was  carried  out  for  a  sufficient  variety  of  conditions /"(0) 
and  ^"(0)  so  that  the  correct  set  of  initial  conditions  yielding 
the  proper  behavior  at  X=  »  (boundary  condition  (12a)) 
could  be  inferred  to  the  desired  degree  of  accuracy. 

Integration  was  performed  according  to  the  following  basic 
scheme:  With  the  value  of  /'"(X)  and  ^'"(X)  given  at  five 
closely  spaced  values  of  X,  fourth-degree  polynomials  may  be 


passed  through  the  two  sets  of  values  of/'"  and  i p"' .  Then, 
if  U  f,  f" ■  V r  and  <p"  are  known  at  the  fifth  point,  the 
polynomial  representations  of/"'  and  may  be  integrated 
to  yield  /  /',  /",  ,  and  at  the  next  (sixth)  point. 

These  quantities  may  then  be  substituted  into  differential 
equations  (11)  to  yield/'"  and  at  the  sixth  point.  In 
this  way,  the  solution  may  be  extended  one  step  at  a  time,  in 
each  step  by  use  of  the  solution  at  the  five  previous  points. 
In  order  to  begin  this  procedure,  the  solution  must  first  be 
found  at  five  points  starting  at  the  origin  and  must  be  subject 
to  boundary  conditions  (12b)  and  (12c)  and  the  tentative 
selection  of /"(())  and  ^"(0). 

This  preliminary  calculation  was  done  in  the  following 
manner:  /'"( 0)  was  calculated  directly  from  equations  (11) 
and  (12)  and  was  used  as  an  initial  estimate  of/'"  at  the  next 
four  points.  Given /(0),/'(0),  and /"(()),  the  values  of/,/', 
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and/"  were  computed  at  the  second  point  by  integrating  a 
fourth-degree  polynomial  passed  through  the  five  values  of 
/'".  In  a  similar  manner,  <p,  p' ,  and  p"  at  the  second  point 
were  found.  Direct  substitution  into  equations  (11)  then 
yields  improved  estimates  for  f"  and  p'"  at  the  second  point 
and  thus  an  improved  polynomial  representation  of  these 
functions  which  may  be  used  to  obtain  values  of  /,  /',  /",  p, 
p',  and  p"  at  the  second  point,  and  so  forth,  until  improved 
values  have  been  obtained  at  the  fifth  point.  This  procedure 
was  repeated  in  an  iterative  manner  until  convergence  was 
obtained  at  each  of  the  five  initial  points. 

All  calculations  were  performed  on  the  IBM  Card  Pro¬ 
grammed  Electronic  Calculator.  Results  are  considered 
correct  to  four  significant  figures. 
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